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We present an overview of a perturbative-kinetic approach to jet propagation, en¬ 
ergy loss, and momentum broadening in a high temperature quark-gluon plasma. 
The leading-order kinetic equations describe the interactions between energetic 
jet-particles and a non-abelian plasma, consisting of on-shell thermal excitations 
and soft gluonic fields. These interactions include 2 2 scatterings, collinear 

bremsstrahlung, and drag and momentum diffusion. We show how the contri¬ 
bution from the soft gluonic fields can be factorized into a set of Wilson line 
correlators on the light cone. We review recent field-theoretical developments, 
rooted in the causal properties of these correlators, which simplify the calculation 
of the appropriate Wilson lines in thermal field theory. With these simplifications 
lattice measurements of transverse momentum broadening have become possible, 
and the kinetic equations describing parton transport have been extended to 
next-to-leading order in the coupling g. 


1. Introduction 

The suppression of highly energetic jets (or jet quenching) is one of the most striking 
findings of the experimental program of heavy-ion collisions.A comprehensive 
review® can be found in this volume, which also contains another contribution 
reviewing a specific approach in greater detail.^ 

In this review, we concentrate on a weakly-coupled kinetic approach describ¬ 
ing the propagation of high momentum jet-like particles through a Quark-Gluon 
Pasma (QGP). A detailed perturbative description of the QGP and jet-quenching 
is available when the temperature is high T ^ Aqcd, and the momentum of the 
jet-particles is much larger than the temperature, T. 
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The kinetic picture of the high temperature QGP, which has emerged through 
a combination of physical intuition, direct diagrammatic analysis, and a gradi¬ 
ent expansion, is characterized by hard particles and random classical non-abelian 
fields.We will describe a Boltzmann equation for the propagation of jet-like 
particles interacting with the hard particles and random fields which comprise this 
idealized weakly coupled plasma. Here the jet-particles have momentum p ^ T, 
while the hard particles have momentum ^ T, and the random classical fields are 
soft with momenta of order gT, p ^ T ^ gT, where g is the strong coupling con¬ 
stant. For simplicity, we will limit the discussion to pure gauge theory. As discussed 
in Section 2, to leading order in the coupling constant there are three processes rel¬ 
evant to the transport of jet-particles in the high temperature plasma: (i) 2 o 2 
elastic scatterings with hard particles, (ii) collinear radiation, and (iii) drag and 
momentum dilfusion driven by the soft classical background. 

The interactions between the jet-particles and the classical non-abelian fields, 
requires a resummation scheme known as the Hard Thermal Loop (HTL) effective 
theory,®’® which is the QCD analog of the Vlasov equations.^® This effective theory 
is necessary to compute the drag coefficients and collinear bremsstrahlung rates 
which describe the propagation of jets. The computational complexity of the Vlaslov 
equations would, at first sight, make any extension beyond leading order in the 
coupling extremely challenging. However, Hard Thermal Loops correlators (and 
statistical correlators more generally) simplify greatly when evaluated at lightlike 
separations.^^ It is precisely such light-like correlators which must be evaluated to 
determine the drag, diffusion, and collinear bremmstrahlung rates of jet-particles 
propagating in plasma. 

These simplifications are a consequence of the following physical picture: Since 
the hard and jet-particles are propagating almost exactly along the light cone, they 
are probing an essentially undisturbed plasma, at least as far as the soft classical 
background is concerned. Informally, we say that the soft classical background 
“can’t keep up” with the hard or jet-particle traversing the plasma. Thus, the soft 
correlations that are probed by the hard and jet-particles are statistical in nature 
rather than dynamical. 

In Sec. 3 we will illustrate the basics of Hard Thermal Loops and then proceed to 
introduce the field-theoretical arguments for the light-like simplifications, originally 
due to Caron-Huot.^® We then use these simplifications at leading order to deter¬ 
mine the transverse momentum diffusion and drag coefficients for high momentum 
probes. 

An important consequence of these lightlike simplifications and their physical 
origin (see above) is that lattice techniques can be used to compute several lightlike 
correlators, which arise in jet-quenching physics. Due to the Minkowskian nature 
of the problem and the large energies at play, a direct ab-initio lattice calculation 
of jet propagation in medium is not feasible. However, as we show in detail, the 
soft contribution to quantities such as q (the jet quenching parameter) becomes 
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Euclidean and three-dimensional, which makes it amenable to lattice measurements. 
In Sec. 4 we illustrate the basic principles of the dimensionally reduced (three- 
dimensional) effective theory, and review the hrst lattice calculation^^ of the soft 
contribution to q. 

Of course, the special properties of HTLs on the lightcone also simplify per¬ 
turbative calculations. Indeed, Caron-Huot employed these results to compute q 
and the transverse scattering kernel to NLO^^ in the strong coupling constant g. 
Building on this and related developments, the photon and low-mass dilepton emis¬ 
sion rates were also computed to NLO.^^d^ These electromagnetic rates at NLO 
correct the strictly collinear approximations made at leading order. Finally, this 
review anticipates an upcoming publication by the current authors where the LO 
Boltzmann equation for jet-transport described above is extended to NLO.^® In 
Sect. 5 we sketch the Boltzmann equation to this order, which includes the drag 
coefficient at NLO and semi-collinear emission rates describing bremsstrhalung at 
wider angles. Many of the parameters of the effective kinetic theory, such as drag 
coefficients, thermal masses, and the transverse collision kernel, can be computed 
using the three-dimensional lattice described above, providing a tantalizing semi- 
perturbative description of energy loss. These results and future directions are 
summarized in Sec. 6. 

2. The kinetic picture at leading order 

In this section we will summarize the kinetics of high temperature weakly coupled 
non-abelian plasmas. As we mentioned in the introduction, this kinetic picture 
of the QGP is characterized by hard particles and random classical non-abelian 
fields.First we will describe the interactions between these constituents quali¬ 
tatively, and subsequently give a more quantitative description. We are particularly 
interested in formulating a Boltzmann equation for the transport of high momentum 
gluons p 3> T traversing the QGP, i.e. a Boltzmann equation for jet-particles inter¬ 
acting with the hard particles and random classical fields which comprise the high 
temperature plasma. The phase space distribution of the jet-particles is notated 
with f{t,x,p) (or fp when clear from context) while the phase space distribution 
of the equilibrium hard particles is notated with n(t,x,p) (or Up). 

The interactions between the particles and the classical helds of the QGP are 
characterized by three processes: 

(1) First, there are 2 o 2 collisions between the hard and the jet-particles. These 
scattering events are strongly localized in space-time. 

(2) Second, the random classical fields in the plasma induce drag and momentum 
diffusion and change the momenta of the jet-particles in characteristic ways. 
The drag and diffusion coefficients are calculated by examining the response of 
the classical fields to the non-equilibrium jet-particle. 

(3) Finally, the momentum diffusion produced by the random classical held induces 
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collinear bremsstrahlung, causing the jet-particles to split collinearly. 


At leading order in the coupling g, the Boltzmann equation for the jet-particles 
takes the form 

/ d d \ 

( ^ + ^ j /(^> p) = C'2«2 [m] + C'diffM + C'coll , (1) 

where the three rates on the right-hand-side reflect these three processes. We have 
anticipated that the 2^2 and diffusion rates depend on a separation scale /r, but 
the dependence on fj, cancels when both processes are included. In the remainder 
of this section we will discuss each of these rates in greater detail. 

The hard particles are approximately on-shell and carry the majority of the 
energy and momentum of the plasma. For this reason these hard constituents are 
the most important for the thermodynamics of the QGP, and the soft fields are 
only important insofar as they influence the kinetics of the jet and hard particles. 
The typical momentum of these excitations is of order the temperature, but the 
virtuality is of order ~ g^T^ ■ Specifically, for a hard particle moving in the positive 
2 direction, the scaling of momentum in lightcone coordinates is®' 

= {p^.P~,P±) , (3) 

= {T,g^T,gT). (4) 


Squaring this momentum shows that the virtuality of these on shell constituents 
is of order ~ g^T^- To estimate the contribution of these on-shell constituents 
to the plasma energy density we use the equilibrium distribution function, np = 
l/(gSp/T _ and write 


e{t, x) = 2d A 



d^p 

(2^)3 ^P’^P- 


( 5 ) 


First, we note that the integral is cut off when the three-momentum of the con¬ 
stituents is soft, p ~ gT. At this point the power counting implicit in Eq. (3) no 
longer applies. Taking the free gas dispersion relation Cp = p, we see that the con¬ 
tribution of the soft modes {p < gT) is of order g^T'^, while the leading contribution 
is of order T^. 

The first correction to the free gas result for thermodynamic properties the 
QGP can be found by determining the correction to the dispersion curve for the 
on-shell particles. Specifically, the dispersion curve (the relation between Ep and p) 
is found by evaluating the one-loop self energy for an approximately on shell gluon. 
The self-energy diagram involves internal lines of momentum 


_ = { q+,q-,q±)-{T,T,T), (6) 

^ We use a notation where p~^ = ^ ^ , p~ = pP — pz, four-vectors are denoted by uppercase 
letters and the metric is the “mostly-plus” one, so that = —2p'^p~ The integration 

measure is 


d^P r dp^dp-d^p± 

(27r)4 “ J 


( 2 ) 
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and are thus highly virtual compared to the on-shell particles. The effect of these 
virtual modes (which are integrated out in any kinetic description) is to change the 
dispersion relation of the on-shell modes 


Ep=p + 



(7) 


where the ellipses indicate higher order corrections. Here we have defined the leading 
order gluon asymptotic mass in pure gauge theory at high temperature^^“^° 




Cpp Ub 
(27r)3 p 


1 

6 


g^CAT"^ ■ 


( 8 ) 


As discussed in Sec. 5, the asymptotic mass can be expressed as a correlator of field 
strengths along the light cone, which is closely related to an analogous correlator 
for q, see Eq. (80). 

These hard modes with p ^ T occasionally collide with the jet partons, transfer¬ 
ring momenta of order the temperature. At leading order these collisional processes 
are 2 -f-)- 2 processes, which can be localized in space-time to within a distance 
of ~ 1/T. A kinetic equation for the phase space distribution function of the jet 
partons including such scattering processes takes the form 


Wt 



fp{t,x) = C' 2 ^ 2 [m] + • • ■, 


(9) 


where the ellipses will be described below. The velocity in this equation is the 
group velocity, Vp = dEp/dp. The collision rate between the hard modes and the 
jet particles takes the form'^ 


~ C ' 2<->2 [ g ] — 


|A^|= 


, , ,, ,, ,, ~(2Try6‘^(P + P'- K - K') 

(l-t-71fc ) (l-t-U;-/) (l-t-TTp/) U-fc/zc'(1-t-Up/) , (10) 


where in the second line we have dropped terms that are exponentially suppressed 
for p ^ T.^ The kinematics of the collisional process with momentum exchange of 
order ~ T is shown in fig. 1. The squared matrix element (summed over initial and 
final state colors and helicities) may be obtained from vacuum perturbation theory: 

|Mp = 16<i.,s‘Ci(3- = -5-(^) . (11) 

The 2 2 scattering rate is divergent in the infrared, and therefore we have 

notated the dependence on the separation scale pt in Eq. (10). There are many 

’^We use a familiar shorthand notation, = f ( 27 r)^ ( 2 ?/)^ ’ neglecting the 

thermal mass when writing Eq. (10). The leading factor of l/(4d^) reflects an average over spins 
and colors of the external gluon and a symmetry factor 1/2 for the sum over final states. 

‘^We consider exp(—p/T) 1. In the remainder of this review we sometimes treat T/P as an 
explicit expansion parameter for illustration purposes. The detailed exposition and the calculations 
in^® are performed without treating T/P as an expansion parameter. 
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Fig. 1. Hard 2 2 collision contributing the collision rate C'2-f^2[M]- hard lines which 

enter or exit the boxed region are included in an effective Boltzmann description. 


ways this divergence can be regulated. At leading and next-to-leading we find 
it convenient^® to simply cutoff the transverse momentum exchange at small q±, 
q± > /r. It is not difficult to extract the logarithmic dependence on ^ for /i ^ 
T. Indeed, let us consider for illustration a leading-log approximation to C2i^2[^A■ 
we expand the distribution function and matrix elements to second order in the 
exchange momentum Q and arrive at a Fokker-Planck equation^ for fp 

■ df 1 ■ ■ (9^ f / T \ 

C 2 « 2 [m] = euvi^J-) + + C> f - j -H /^-independent, (12) 


In writing this equation we have dropped terms suppressed by Tjp. Here t) is a unit 
vector in the direction of p, and the diffusion tensor q^y {p) controls the longitudinal 
and transverse momentum diffusion, 

Quvip) = qL,uv{p)v'‘v^ + ]^quv{pW^ - v^v^) ■ (13) 


The values of these coefficients are found from the expansion of Eq. (10), and for 
pure gauge are at leading log 


Quvip) = 9 ^CaT^ log , (14) 

qL,uv{p) = 9 ^CaT'^ log . (15) 


Here the Debye mass is given by the integral over distribution functions 


ml, = 2g^CA 


d^p np(l + np) _ 1 2 ,^ ^2 

(27r)3 T 3^ ^ 


(16) 


and the asymptotic mass is given by a similar integral in Eq. (8). At this point 
the interpretation of these thermodynamic integrals as the Debye and asymptotic 
masses is premature. This interpretation will be clear from Sec. 3, which explains 
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physically why the transverse diffusion coefficient involves the Debye mass, while 
the longitudinal diffusion involves the asymptotic mass. 

The fluctuation-dissipation theorem relates the coefficient of qL,uv{fJ-) to 

qL,uv{fJ') = + O . (17) 

This Einstein relation relation is necessary in order for fp to reach its equilibrium 
form under the Fokker-Planck evolution in Eq. (12). Since the equilibrium 

distribution is a stationary solution to the full 2 <->■ 2 Boltzmann equation (Eq. (10)) 
this relation between the drag and the longitudinal diffusion coefficient automati¬ 
cally arises from a direct expansion of the Boltzmann equation. 

As we will see, the cutoff dependence in the 2^2 contribution to the Boltzmann 
equation cancels when the interactions with the soft background fields are included. 
The description of these soft modes is our next task. As discussed above, the 
virtuality of the hard modes is of order ^ (gT)^ with the momentum scaling 
given in Eq. (3). The soft gluonic excitations, on the other hand, have a momentum 
scaling 


P>^ = {gT,gT,gT), (18) 

but also have virtuality of order ~ (gT)^- Since these modes are soft, they are 
highly occupied due to the Bose-Einstein distribution function 

J ^ J ■ ( 19 ) 

Thus, at leading order these soft modes can be treated as a classical gluonic field 
interacting with the hard on-shell modes. 

It is important to emphasize that the energy in these soft modes constitutes a 
small fraction of the total energy density. Nevertheless, they significantly influence 
the kinetics of the on-shell particles. Specifically, these soft modes collide frequently 
with the hard and jet particles and exchange soft momenta of order gT. It seems 
intuitive that the effect of these soft scatterings can be incorporated into a Fokker- 
Planck equation describing the drag and momentum diffusion of the jet (and hard) 
particles 

+ Vp ■ fp = C'diff [m] + ■.. , (20) 

The collision (or Fokker-Planck) operator can be written as 

= + ( 21 ) 


where the momentum diffusion parallel and perpendicular to the direction of motion 
p are quantified by and the fluctuation-dissipation theorem is satisfied as 
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described above'^. Terms of higher order in T/p are discussed in detail elsewhered® 
We have emphasized that the parameters of the Fokker-Planck evolution depend 
on the scale p separating the scattering from soft modes, and the scattering from 
hard 2 e-)- 2 collisions. As the separation scale is changed, the parameters of the 
Fokker-Planck equation change and the drag and diffusion rates change accordingly. 
Indeed, Sec. 3 shows that at leading order the drag and diffusion coefficients are 

q{p) =g^CAT^\og f —) , (22a) 

27r \ mo ) 

qUp) =g^CAT'^ log (. (22b) 

27r \m^J 

Thus, a change in these coefficients due to a change in p is compensated by a 
corresponding change in the hard 2 e-)- 2 collision operator in Eq. (12). 

As is well known the Fokker-Planck evolution specified by Eq. (21) is equivalent 
to a Langevin process,where each particle in the phase space follows an equation 
of motion specified by a drag and a stochastic force whose variance is specified 
by the diffusion coefficients 


^ = -e{py+e, (23) 

The Langevin description is valid on timescales which are long compared to the 
underlying correlation times between the microscopic forces in the medium. For 
a classical particle satisfying Newton’s Law, dp'‘/dt = J"*, one would adjust the 
diffusion coefficient so that the integrated squared variance between the microscopic 
and stochastic forces agree 


q^^T= JdtJ dt' (e(t)e (t')) = J • (24) 


Here T is the total time, and the forces are to be evaluated along the trajectory 
of the particle, i.e. for a particle traveling at the speed of light along the z-axis 
the coordinates are = {x^,x~,x±) = (t, 0,Oj_). For a hard particle in repre¬ 
sentation R interacting with a soft classical background QCD field, the forces are 
gT^F^f^{x~^)Vfj_ and are dressed with Wilson lines following the trajectory of the 
particle along the light cone from past infinity'’ 


R'(x^) = U]^{x''~, -oo) gF’-f^{x^)Vfj, Ur{x''~, -oo) . (25) 

^In the remainder of this review q describes only the contribution of soft gT to transverse momen¬ 
tum broadening, and perhaps should be written as ^soft- most of the literature on jet energy 
loss (see for example q refers to the sum of the soft and hard contribution, with a logarithmic 
sensitivity on the jet energy p. 

®Here and below we use a matrix notation is a lightlike four vector = (Ij'u) 

in the direction of motion of the particle, which is conventionally taken to be along the 2 axis. 
Thus, = F^~ and note that F^~ = F'^ . More explicitly, the Wilson lines are f/(ai+,0) = 

Pe~'^ ^0'^ dx+A- _ 
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The appropriate formulas for q and ql are then 


1 

g(^) =— J dx+ Tr(j'j_(x+)J'j_(0)) , 

/ dx+ Tr{F^{x+)r\Q)) , 


J — oo 


— oo 


(26) 


(27) 


where the two transverse directions are summed over in Eq. (26), and the trace 
averages over the colors of the incoming particle. For a classical background field 
these correlators can also be written as^ 



2 /^ roo 

9L(M)=^y_ dx+{v^Fi:+{x+)Ut{x+,0)vpFP+{0)) , (30) 


where U'^{x~^,0) is the adjoint Wilson line. 

The evaluation of the transport parameters qL{^J-) and g(/r) using these correla¬ 
tors involves understanding how the hard particles interact with the classical gluon 
fields at leading and next-to-leading order. This had been carefully examined by 
Blaizot and lancu who worked in a background field gauge and systematically em¬ 
ployed a gradient expansion to determine the appropriate kinetic equations which 
are the non-abelian generalization of the Vlaslov equations.We will write the 
non-abelian Vlaslov equations in the next section and exhibit the appropriate HTL 
diagrammatic rules for evaluating these soft classical gluon correlators. Then in 
Sec. 3 we will show how these rules can be used to arrive at Eq. (22). 

The random classical gluonic fields cause the hard and jet particles to diffuse in 
momentum. This random walk induces collinear bremsstrahlung in medium. We 
will follow the AMY formalism for computing the collinear bremsstrahlung rate,^^ 
and refer to ref.^® which compares the AMY formalism to other approaches, which 
are described in the reviews.®’® Collinear bremsstrahlung can be included into the 
Boltzmann equation and takes the form of a local rate Ccoii for a hard (approxi¬ 
mately massless) gluon to branch into two hard gluons moving in approximately 
the the same direction as their parent^®^®® 


(^t T ' ^x') fp - Ccoll -j- . . . . 


(31) 


To describe the gluon bremsstrahlung rate we refer to fig. 2 which schematically il¬ 
lustrates the collinear radiation induced by the soft scatterings. The kinematics are 
shown for a gluon with momentum p = |p| in the longitudinal direction to branch 
into two gluons moving in approximately the same direction with longitudinal mo¬ 
mentum w and p — uj. Although the bremsstrahlung rate enters the Boltzmann 

^ To see this, take dx'^ = Ax'^ small. Then, evaluate the following trace to first order in Aa;+ to 
see the structure of the adjoint Wilson line emerge 


Tr[J'"^(x+).F^(0)] ~ Vf,F^t^{x+)vpF^'’{0)T[[{l + igAx+A-)T‘^ {1 - igAx+A-)T^] . (28) 
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Fig. 2. Schematic Feynman diagram contributing to the leading order collinear bremsstrahlung 
rate. Hard gluon lines are labeled by their three momentum {pz,p±). The interactions with the 
random classical background bath are illustrated by the gluon lines with crosses. Only hard lines 
which enter or exit the boxed region are included in an effective Boltzmann description. 


equation as a local rate, it must be understood that the emission process can only 
be localized to within a time scale set by the formation time of the radiation. The 
inverse formation time will be defined as the energy difference between the initial 
and final states 


(’^form) — 5E(h^p^uj) — T Ep—^^ Ep . 
Using the dispersion relation for the hard particles this reads 


SE{h,p, uj) ~ 


'^’ooui "^oop-bj 


OO p 


2poj{p — oj) 2uj 2{p — oj) 2p ’ 


(32) 


(33) 


where ^ is the asymptotic mass of the particle with momentum p, as summarized 
in Eq. ( 8 ). We have further defined 


h = pq± . 


(34) 


As seen from the figure and described below, h/p is a transverse momentum vector 
which is conjugate to the (transverse) coordinate separation xj^ between the initial 
and final states. 

The bremsstrahlung rate Ccoii is determined by the rate of transverse momen¬ 
tum kicks (of magnitude q±) which a hard particle experiences traversing the soft 
classical fields: 

Cfl(<?J-) = lim , ( 35 ) 

P-^oo d‘^q± 

Here p is the momentum of the hard particle, which is large {p —^ 00 ) relative to 
the the typical momentum, ~ gT, of the background fields. The collision kernel 
Cr can be expressed as a Wilson loop in the (a;“'",x_L) plane evaluated in the clas¬ 
sical background,as sketched in Fig. 3. To motivate the appropriate Wilson 
loop we note that the average squared momentum transfer per unit time (i.e. q) is 
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determined by Cji{q±) 

9(Ai) = J . (36) 

As described above q can be expressed as a correlation of field strengths (see 
Eq. (26)). Numerically at least, this correlation function should be understood 
as the limit of a Wilson loop (see fig. 3) as the transverse separation approaches 
zero and the length approaches infinity 

lim lim (WA(a;_L,L)) = exp (-^q(^)a;^L ) (37) 

Si—>0 L—s-oo V / 

It is not difficult to show that this Wilson loop expression agrees with Eq. (26) 

s+ (0.0) (L.O) 


-^ 

-> 

V 

<- 


(O.ii) IL,X±) 

Fig. 3. Wilson loop leading to g(/r). Here the side rails, Pe“*Io (s+.o) 

^ are separated by transverse coordinates Xi and extend in a:+ = (t + 2)/2 

from a:+ = 0 up to a;+ = L. x~ is zero in this figure. The Wilson loop is traced over colors, 
^(^Tr..., as indicated by the gap at the top left corner. 


in the limit of small separation. Although it is not entirely obvious, it is certainly 
not surprising that the Wilson loop is sufficient to completely determine Cii{q±). 
Specihcally, the impact parameter representation of the collision kernel 

Cr{x^) = I ^ (1 - CR{q ^), (38) 

is determined by the asymptotic behavior of the same Wilson loop at large L 

lim {Wr{x±,L)) = exp{-CR{x±)L) . (39) 

L—¥00 

Expressions of this form have been used to compute g(/r) at leading and next-to- 
leading order. 

We have described the collision kernel in detail, since it is a key parameter in 
determining the collinear bremsstrahlung rate. A derivation of this rate is beyond 
the scope of this review, and here we will simply quote the relevant formulas and 
describe the relevant physics. The radiation rate for a gluon of momentum p to 
radiate a gluon of longitudinal momentum w in approximately the same direction 
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IS' 


5,31 


dT{p,u}) 

du! 


coll 


9^C!a , , , 1 + + (1 — 

j^(l + ,.(■.))(! + n(p - 


d^h 


2h • ReF(h,p, w), 


(40) 


where x = ui/p is the longitudinal momentum fraction of one of the outgoing gluons. 
The leading factor (1 + + (1 — x)^)/(x^(l — x)^) records the collinear emission 

probability accompanying a scattering process. Similarly F(h, p, oj) records the time 
evolution of the current over the formation time. In momentum space this time 
evolution results in an integral equation for F(h,p, w) which resums the influence 
of multiple soft scatterings® on the emission process®’^^ 

2h = zSE(h,p,uj)F(h) + J |^^^M|[F(h)-F(h-a;k^)] 

+ [F(h) - F(h + pkx)] + [F(h) - F(h -{p- a;)kx)]|. (41) 


Finally, the full contribution to the Boltzmann equation takes the form of a gain 
term (reflecting the radiation from particles with momentum p + uj) and a loss term 
(reflecting the radiation from the particle with longitudinal momentum 


l^coll — 


duj fp-\-uj 


dr(p + w, Lu) 


du; 


- 9{p - 2u;)fj 


dr{p,u;) 


coll 


duj 


(42) 


coll 


The region with w < 0 represents the collinear merging of two gluons to make a final 
state gluon. We also remark that Eq. (42) is finite in the infrared, as the power-law 
divergences for a; —>■ 0 and w —>■ p cancel between the gain and loss terms. 

In practice, the integral equation for the bremsstrahlung rate (Eq. (41)) is most 
easily solved by working in coordinate space, which also illustrates the role of the 
coordinate space Wilson loop in shown in fig. 3. To facilitate this rewriting, we recall 
h/p is conjugate to a;_L, which we relabel as b. Dividing Eq. (41) by p and Fourier 
transforming with respect to h/p we arrive at a coordinate space representation of 
the integral equation: 


- 2iVS'^{h) = 


- X + x^) - \7l) f{h) 


2x (1 — x) 

fCAixh) , CA{h) , C^((l-x)b) 






/(b), 


where we have defined 


/(b) = 


drh , 

rC p F(h) 


(27r)2p3 


(43) 


(44) 


^Their effect on the emission process is usually called the Landau-Pomeranchuk-Migdal (LPM) 
effect. 
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and note that the integral entering the rate (40) becomes 

^ 

This differential equation (Eq. (43)) can be solved to determine the bremsstrahlung 
rate. 

The coordinate representation of the integral equation (Eq. (43)) naturally in¬ 
volves the Fourier-transformed collision kernel Ci^(b) stemming from the Wilson 
line discussed above (see fig. 3 and Eq. (39)). Let us remark that the second line 
in Eq. (43) contains a sum of three two-body contributions, with the transverse 
distance b weighted by the ratio of two of the three energies p, uj and p — oj. Indeed, 
it is easy to show^^’^®“^® that, at leading order, this is the only possibility. Beyond 
leading order, one could expect a more complicated three-pole structure, with three 
Wilson lines on the light cone, properly linked by transverse lines at their endpoints. 
However, at least up to NLO the symmetries of the system yield the simple planar 
dependence in the transverse plane reflected in Eq. (43).^^ 

To summarize, we have identified three processes at leading order which are 
relevant to the transport of high momentum jet particles in a weakly coupled QGP. 
The first process is straightforward 2 o 2 collisions. This rate depends on the sepa¬ 
ration scale /i between the ^ T scale, characterized by hard particles, and the ^ gT 
scale, characterized by classical fields. The second process is momentum diffusion, 
which describes the interactions between the classical fields and the jet-particles. 
This processes is described by a Fokker-Planck equation, whose parameters e(/i) 
and q{p) also depend on the separation scale. However, the sum of the 2 <->■ 2 
collisions and the Fokker-Planck evolution is independent of p. Finally, the third 
process is collinear bremsstrahlung, which is independent of the separation scale 
to this order. The parameters entering the collinear bremsstrahlung rate are the 
asymptotic masses and the scattering kernel C{q±). Both of these quantities 
can be expressed as (closely related) correlators of field strengths along the light 
cone, Eqs. (29),(39), and(80). 

Anticipating Sec. 5, at next-to-leading order we will find corrections to the pa¬ 
rameters of the Fokker-Planck equation, e{p) and q{p), and corrections to the pa¬ 
rameters of the collinear bremsstrahlung rate, and C{q±). As we will see the 
NLO corrections to e(p) depend logarithmically on the separation scale p. How¬ 
ever, at NLO a correction to the collinear bremsstrahlung rate must also be included, 
which accounts for not-so-collinear emissions. This semi-collinear emission rate also 
depends on the separation scale p, and the complete NLO rate is independent of p. 
Thus, the NLO correction to the collision term takes the schematic form 

5C = ^Ccoll + -I- dCsemi-coll M ■ (46) 

Before we describe these developments in detail, we will show how the compu¬ 
tation of lightlike correlators in thermal field theory (such as e{p) and q{p)) can be 
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dramatically simplified using the causal properties of these correlators. This is the 
task of the next section. 

3. Hard Thermal Loops on the lightcone: a modern perspective 

In the previous section we described a kinetic equation for jet-like particles based 
on 2 o 2 collisions, diffusion and drag, and collinear bremsstrahlung. The diffu¬ 
sion and drag coefficients are determined by classical statistical correlators of the 
soft background held on the light cone, Eq. (26) and Eq. (27). These correlation 
functions are to be evaluated using the Hard Thermal Loop (HTL) effective theory, 
which describes the response of the coupled Maxwell-Boltzmann system. Similarly, 
the collision kernel C[i{q±) is also determined through Eq. (39) by the interaction 
between the classical particles and the the background held on the light cone. 

Hard Thermal Loops involve a complex set of medium-modihed propagators 
and vertices. Already at leading order the evaluation of the necessary correlators is 
seemingly complex. However, as we will describe in this section the evaluation of 
lightcone HTL correlators turns out to be surprisingly simple. For certain operators, 
such as those associated with transverse momentum diffusion Eq. (26), and soft 
corrections to the asymptotic mass, a Euclidean formalism can be used.^^ On the 
other hand, other operators such as those associated with longitudinal dihusion and 
drag are surprisingly insensitive to the ^ gT sector, and lead to contact terms which 
correct the collinear emission. 

These simplihcations are a consequence of the following physical picture. Since 
the hard and jet-particles are propagating almost exactly along the light cone, they 
are probing an essentially undisturbed plasma, at least as far as the soft classical 
background is concerned. Thus, the soft correlations that are probed by are the hard 
and jet-particles are statistical rather than dynamical, i.e. these are correlations that 
arise from the partition function, rather than a retarded response averaged over the 
correlated initial conditions determined by the partition function. This heuristic 
explanation underlies the dimensionally-reduced lattice calculation of q which is 
described in Sec. 4. 

At a technical level, these simplifications were first realized using light cone 
causality of retarded propagators,^^ which says that 

^ (47) 

is an analytic function of p'^ in the upper half-plane at fixed p~ and p±. This is 
because the retarded response function is only non-zero in the forward light cone 
2x~^x~ > x^. Thus the integral in Eq. (47) has support only for x~ > 0, and the 
Fourier integral provides an analytic continuation in the upper halfp"*" plane, due to 
the decreasing exponential “ . To illustrate these simplihcations in a concrete 
context, we will evaluate the drag and transverse momentum diffusion coefficients, 
Eq. (26) and Eq. (27), at leading order. To do this we will hrst briefly review Hard 
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Thermal Loops in Sec. 3.1. Then we will use the HTL theory to evaluate Eq. (26) 
and Eq. (27) in Sec. 3.2 and Sec. 3.3, using light cone causality to simplify the 
calculation. 

3.1. Hard Thermal Loops and Light Cone Correlation Functions 

To briefly describe the Hard Thermal Loops consider a scalar field (j)a transform¬ 
ing in the representation R (where a is a color index). There is a single particle 
density matrix Ni,a{t,p, x) transforming in the R x R representation which records 
the statistics of the color orientations of a quasi-particle excitation. Heuristically,^® 
Nba{t,p,x) is a quasi-local expectation value of (^a\ p{t, x)ab,p{t, x)') at the space 
time point {t,x). More formally, Nba{t,p,x) can be defined as the Wigner trans¬ 
form of the gauge-covariant Wightman function. The gluon singlet distribution 
function (denoted Up above) is proportional to the color trace of a single particle 
density matrix. 

Tip = ^Tr [N(t, X, p)] (48) 

UR 

while the color current associated with the gluon distribution is** 

J^^^{t,x) = 2g I I^Tr [T^N{t,x,p)] (49) 

where = (l,p). The single particle density matrices Nba obey a Vlaslov equation 

{Dt+v- D^)Np +^gS^{E + vx B)z, = 0, (50) 

where Dp = dp — ig{Ap, •] denotes the appropriated gauge covariant derivative for 
the R X R representation, and the curly brackets denote an anti-commutator. The 
particles source the classical chromo-electric and magnetic fields 

-DpF>^''{t,x) = r{t,x). ( 51 ) 

The Hard Thermal Loop retarded Green function is derived from this Maxwell- 
Boltzmann system of equations, Eq. (50) and Eq. (51), by expanding around the 
equilibrium state, Up — l/(e^p/^ — 1) 

Nba =np6ba + 5Nba , (52) 

Ap =0 -l- 5Ap , (53) 

and solving for the response of the classical field to an external current, A>^(uj, q) = 
iC^{u},q)J^^ext{i^,q)- We will not go through this set of steps but refer to the 


*^The leading factor of 2 is the spin of the gluon, and for a scalar field this factor would not appear. 
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literature.The resulting retarded Green functions in the Coulomb gauge are 


q^ + TlL{q°/q) 




Gl{u,q) = 


- q'q^) 


_(,^0)2+q2+nT(g79) 

where = — 1 and the self-energies are 

a: -I- 1 


q^—uj+ie 


nL(a;) =mjj ( 1 - - log 


X —1 


m 


T\t{x)=^{x^- 


{x^ — l)a: /a; -I- 1 


log 


a; — 1 


(54a) 

(54b) 

(55a) 

(55b) 


with X = q^/q. The spectral density is the difference between the retarded and 
advanced green functions 


pf^''{u})=G'^iu},q)-G'^{ui,q), (56) 

where GA{uJ,q) is found by setting = oj — ie in Eq. (54). The HTL spectral 
density and its properties are discussed in the references. 

As we described above, q and (Jl are determined by the propagation of the 
hard and jet particles in a random classical background field. In equilibrium, 
the fluctuation-dissipation theorem relates the symmetrized two point function 
Grrixi, q) to the imaginary part of the retarded response' 

G(l,^(w,g)^i({A^(a;,9 ),A''})= (^nB(w) + i)p^'^(a;,g). (57) 

For a classical background field the symmetrization of the field operators is 
irrelevant, and the symmetrized two point function simply records the two-point 
statistics of the classical field. Approximating Eq. (57) for a; ^ T yields 

G(f,^(a;, q) = (^'^(a;, q)A'') = Jw, q) . (58) 

UJ 

We will use these statistics of the classical soft modes when evaluating the longitu¬ 
dinal drag and diffusion coefficient in the next section. 

In addition to these rather complex propagators, the Maxwell-Boltzmann system 
leads to additional HTL vertices which describe how the soft classical background 
interacts with itself. The resulting Feynman rules describing this response are rather 
complex as well. To date, this elaborate set of Feynman vertices has been used to 
compute the dilepton emission rate at small invariant mass^^ and the diffusion coef¬ 
ficient of a heavy quark at NLO.^^ As discussed above and in the next sections, these 
intricate HTL rules simplify dramatically when evaluating the lightcone correlators 
involved in energy loss. Indeed, the NLO computations of energy loss described in 
Sec. 5 would not have been possible without these simplifications. 

^Our notation Grr follows the one of the so-called r/a formalism of Thermal Field Theory. We 
refer to this reference^^ for a comprehensive set of Feynman rules for the HTL theory in that 
formalism. 
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3.2. Transverse diffusion and Euclidean operators 

In this section we will compute q and C{q±) at LO using Eqs. (26), (35) and (39) 
as our starting point. Recall that these Wilson line correlators describe the prop¬ 
agation of a hard or jet particle traversing the classical field provided by the soft 
modes. The statistics of these classical fluctuations are given by Eqs. (57) and (58). 
As we will see, light cone causality Eq. (47) dramatically simplifies the calculation, 
leading to a Euclidean formulation. 

Eor C, applying the definition (39), one finds (in any non-singular gauge) that the 
diagrams in Eig. 4 contribute. Clearly, the diagram on the right and its symmetric 



Fig. 4. Diagrams for the leading-order evaluation of C{x±) in a non-singular gauge. The plain 
lines represent the Wilson lines along the light-cone, at transverse separation x±. The diagram 
with the self-energy on the upper Wilson line is not shown. 

give rise to the a;_L-independent part in Eq. (38), whereas the one on the left gives 
rise to the a;_L-dependent exponential, so that, at leading order, for an adjoint source 
(a gluon) one has 

C{qj_) = g^CA J ^^27rS{q^ - q.)G--{Q), (59) 

where G~~{Q) = v^v^GI^^ (Q) is given for soft momenta by Eq. (58). Since = q^ 
enforces space-like kinematics, it would then seem that one is left with an integral of 
the non-trivial HTL spectral function Eq. (56) over the entire Landau cut, < q. 
Such integral poses no particular challenge numerically,^^ but any extension beyond 
leading order, with loops composed of HTL propagators and vertices generated by 
the rules discussed in the previous section 3.1, would be extremely challenging 
from the computational viewpoint. We refer to the computation to NLO^^ of k, 
the analogue of q for a very massive quark, for an example of such an intricate 
calculation. In the present case, however, there is a considerable simplification. As 
we have remarked in Sec. 2, C{q±) can be obtained from a Wilson loop supported 
on a light front (a null plane), as show in Fig. 3. The special analiticity properties 
of objects of this kind are the key to their simplification. Indeed, Alves, Das and 
Perez,and Weldon,introduced the technique of light-front quantization into 
Euclidean Thermal Field Theory, extending the usual imaginary time formalism 














February 13, 2015 1:25 


World Scientific Review Volume - 9.75in x 6.5in 


eloss 


18 


J. Ghiglieri and D. Teaney 


from the space-like domain to light fronts. Subsequently, and independently, Caron- 
Huot^^ introduced a similar quantization technique and showed how it can reduce 
the calculation of (the soft contribution to) this Wilson loop to a much simpler 
Euclidean calculation in dimensionally-reduced EQCD. 

The first step in this endeavor consists in showing that the calculation of n- 
point correlation functions, where all fields lie on a spacelike hypersurface, can be 
carried out by Euclidean techniques.The null correlators we need can also 
be computed provided that they are free of collinear singularities, which they are 
up to NLO included.The proofs presented in Refs.^^’"^®’^® are rather complex and 
technical. Here we will present a much simpler derivation, also due to Caron-HuoE, 
already collected in^^ , which we follow closely. It is an illustration of the two-point 
function case, which, as one can see from Eq. (59), will nicely yield the LO result 
for C{q±). 

Consider the symmetrized correlator of some operator (such as G in our case), 
Grr{x°,x) with \x^\ > |x°|. Since the separation is spacelike, operators commute, 
and therefore Grr equals G^, or G^, as shown in Eq. (58). We are interested in 
the lightlike limit There are additional subtleties in taking this limit, but 

these are irrelevant for the soft and leading order contributions to such correlators. 
Let us write Grr in terms of its Fourier representation, i.e. the inverse of Eq. (47) 

Grr{x°,x) = {2Tr)-‘^ J doj J V G^,.(w, p^p_L ), (61) 


and use the exact relation between the spectral function and the symmetric propa¬ 
gator, Eq. (57). We then shift one integration variable to — {x^jx^)uj: 



dujdp^d^px 

(27r)4 



^ Pi'Xl.p'' +Uj{x'^/x^),px ), 


(62) 

so that the Fourier exponent is now frequency-independent and we can perform the 
LO integration by contour methods. The spectral weight is written as 


p{io,p) = Gr{oj,p) - Ga{lo,p) , (63) 

and the integral over Gr{lo,p) is performed by closing the contour in the uppper 
half plane, while the integration of Ga{i^,p) is performed by closing the contour in 
the lower half plane. Since \x'^/x^\ < 1, Gr(uj,p^ + uj(x° /x^),p±) is an analytic 
function of lo in the upper complex lo plane. This is a consequence of (slightly 
generalized) light-cone causality'^ as described in and after Eq. (47). The advanced 
function is similarly an analytic function in the lower complex uj plane. Therefore 
the only singularities encountered in continuing the frequency integration are those 
in the statistical function (nB(w) -I- 1/2), which has poles at = 2'KinT with 

JS. Caron-Huot, oral presentation at the Institute for Nuclear Theory (Seattle), 29 March 2012 
^ This is precisely a statement of lightcone causality as described in Eq. (47) when jx^ = 1. In 
this case, = —p~ and a; parametrizes the analytic continuation in p"*" at fixed p“ and p±. 
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n = (... — 1,0,1,...) and residues equal to T. Closing the contour around these 
poles, and renaming to , we find' 


Grr{x^, x) = T 




d^p 

(27r)3 


G^PGE[u:r.,P^+iUJn — ,P± 


oJn = innT , (64) 


where we have recalled that the retarded green function is determined by an analytic 
continuation of the Euclidean function, G/j(a;,p) = —iGE{<-OE = —*(w + ie),p). 

When we need to compute the soft gT contribution to such a correlator, one may 
drop the nonzero Matsubara frequency contributions; that is, we keep only the n = 0 
term in the sum. This term, GE{<-OmP^ + /x^)iP±) —t Ge(^,P^,P i.), is the 

Euclidean correlation function of the 3-dimensional dimensionally reduced EQCD 
theory, which we shall describe in more detail in Sec. 4. Thus, the soft contribution 
to Euclidean operators at spacelike and lightlike separations is time-independent. 

Let us go back to the computation of C(a::j_). As all points in that Wilson loop are 
either at a space-like or light-like separation, the above introduced EQCD reduction 
is applicable. Hence, at leading order, Eq. (59) becomes 


C{qE)=9^GAT [ ^27T6{q,)G-E-{0,q,,q±) = O^GaT (^ ^ ^ A , (65) 

J [^TTj \g_L q^+rrijjj 

The result is easily understood as the difference between the massless, transverse 
and massive, longitudinal EQCD propagators, which also arise as the q° = q^ = 0 
limit of Eq. (54). For (the soft contribution to) q we then have the well-known 
logarithmic dependence on the cut-off, which we anticipated in Eq. (22) in Sec. 2, 
i.e. 


dip) = 9 ^CaT j 






m 


D 


p'^GaTw? 


D 


{2t:Y q\{q\+rri]j) 2 tt m_D ' 

Furthermore, we remark that Eq. (65) was already obtained^® using a sum rule. 
Caron-Huot has shown'^ that the two approaches are equivalent, and that the sum 
rule is possible because of the aforementioned analyticity of retarded and advanced 
functions. These very same properties will be of great importance for the evaluation 
of non-Euclidean operators. 

Finally, let us consider the effect of the non-zero Matsubara modes we have 
neglected in Eq. (65). If g_L ^ gT is kept soft, the subtleties in putting jx^ = 1 
will need to be reexamined when including the non-zero Matsubara modes in the 
sum. The evaluation of these contributions is related to the collinear singularities 
that should show up at NNLO. It would be interesting to clarify these collinear 
singularities with direct calculations. If one considers shorter transverse separations 
(with q± ~ T) then all Matsubara modes will contribute straightforwardly, and the 
Euclidean light-front computation in Eq. (64) may be more cumbersome than a 
direct computation in the real-time formalism. The leading-order contribution to q 
from the scale T has been computed^® using real time methods. 

* The pole at n = 0 is an artifact of the separation of p, which vanishes for a; = 0, into Gji and Ga- 
The individual poles there can then be dealt with in a principal value prescription, for instance. 
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3.3. Longitudinal diffusion and non-Euclidean operators 

As we mentioned at the beginning of Sec. 3, not all lightcone or light-front supported 
operators admit a three-dimensional, Euclidean description for the soft modes. A 
prime example is the longitudinal momentum diffusion coefficent q^, as given by 
Eq. (27). At leading order it is given by the diagram shown in Eig. 5. In any 



Fig. 5. The leading-order soft contribution to q^. The two dots are the two field strengths and 
the double line is the adjoint Wilson line connecting them. The curly line is a soft HTL gluon. 

non-singular gauge it reads 

= g^CA dx+ J (67) 

where again G~~{Q) is given by Eq. (58). The x~^ integration sets q~ to zero. 
We clearly see that, although originating from a lightcone operator, q'^ cannot be 
evaluated in EQCD: indeed, the zero-mode contribution exactly vanishes when the 
previous techniques are applied. 

We can however evaluate Eq. (67) by employing sum rules that are rooted in the 
same analyticity properties that were used in the derivation of Eq. (64). In detail, 
we plug the result of Eq. (58) in Eq. (67). Up to 0{g^) corrections™ we then have 

qL=g‘^CA j ^££^Tq+{G-,^-{q+,qj_)-G^-{q+,q_L)). ( 68 ) 

This too would be a simple enough numerical integral® over the HTL spectral 
function in the Landau cut, of difficult extension to higher orders. However, as 
we have previously remarked, retarded (advanced) two-point functions are analytic 
in the upper (lower) half-plane in any timelike or light-like variable. We can thus 
deform the integration contours^® away from the real axis ontoCfl (|9'^| = /r’*' S> gT, 
Im( 7 + > 0) and Ca (|9'''| = ^ gT, Imq+ < 0), as depicted in Eig. 6.*^ /r+ is a 

"^When expanding the statistical factor in the soft region in Eq. (58), one has nB(w) + 1/2 = 
T/uj{l + Oig^)}. 

*^The longitudinal and transverse contributions to (Q) contain poles at g"*" = q~{2 ± iq± 

which, being on both sides of the complex plane, appear to violate analyticity. However 
their residue cancels in the sum of longitudinal and transverse components. As observed in^^ , 
they are artifacts of the decomposition into Lorentz-variant longitudinal and transverse modes and 
their contribution has to vanish in all gauge-invariant quantities. 













February 13, 2015 


1:25 


World Scientific Review Volume - 9.75in x 6.5in 


eloss 


Parton energy loss and pT broadening at NLO in high temperature QCD 21 



Fig. 6. Integration contour in the complex q'^ integration, and the deformation we use to render 
q'^ ^ gT. Gfi runs above the real axis and Ga below. 


UV cutoff on 5 +, to make sure that we consider only the soft region. As we shall 
see, it will have no effect. 

Along the arcs q'^ is then large, albeit complex. The longitudinal and transverse 
propagators simplify greatly, i.e. 

= G^{q+,q_L) + ^G|J(g+,gj_) (69) 

{q+y ^ {q+yql+ml^^^{{q+y 

where we have used the Coulomb gauge propagators given in Eq. (54) and rn^ = 

/2 is the LO gluon asymptotic thermal mass as given in Eq. (8) . What this equa¬ 
tion is telling us is that, for large (complex) g"*", the longitudinal part is unchanged 
with respect to its unresummed version, whereas the transverse part acquires the 
asymptotic mass, which is typical of collinear, transverse excitations. 

Upon plugging Eq. (70) in Eq. (68), the g+ integration along the arcs is trivial 
and the end result is^® 







9^GaT 


(71) 


(27r)2 ql + 27r 

where contributions smaller than l/(g+)^ in Eq. (70) are not needed, as they would 
only give rise to power-law terms in the cutoff on q'^, which would then cancel 
against contributions from larger scales. The transverse integration has been regu¬ 
lated in the UV with fi ^ gT, as in Eq. (66). Again, this result had been anticipated 
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in Eq. (22). As we commented after Eq. (16), this calculation, together with the 
one in the previous subsection, clearly shows the physical reasons causing the Debye 
mass rriD to appear in q{fj,) and the asymptotic mass rrioo to appear in gL(p.). 

The sum rule we have just illustrated has been derived in this fashion in^® . It 
can also be seen as the bosonic equivalent of the one presented in the calculation of 
the thermal photon rate to NLO^^ and of the thermal ultrarelativistic right-handed 
neutrino rate at LO.®° We also remark that the same result (71) was previously ob¬ 
tained in a different way®^ for the drag coefficient e, which is related by the Einstein 
relation Eq. (17). As it is shown there,once the difference in regularization be¬ 
tween q_\_ < /i_L and q < is taken into account, Eq. (71) agrees with the numerical 
results of Braaten and Thoma"^® for r; —>■ 1. 

In summary, in this Section we have introduced Hard Thermal Loops and their 
simplification on the light cone. In Sec. 3.2 we have analyzed the subset of oper¬ 
ators, including transverse momentum diffusion, that admit a three-dimensional. 
Euclidean description in the soft region, whereas in the current subsection we have 
analyzed those that do not, showing how they still considerably simplify. In the 
following Section we will show how Euclidean operators can be calculated on the 
lattice, whereas in Sec. 5 we will use the results of this Section to generalize the 
perturbative kinetic approach to NLO. 


4. Non-perturbative determinations of C{x±) and q 

As we mentioned previously, the field-theoretical developments introduced by 
Caron-Huot^® make it possible to map certain operators supported on light fronts 
and light cones to the Euclidean domain, where the soft contribution is dominated 
by the dimensionally reduced zero modes of 3D EQCD. Recently, this has been 
exploited to compute the soft contribution to C(xi_) and q on the lattice. Laine 
and Rothkopf®® have performed a first study of the Wilson loop defined in Fig. 3 in 
classical gauge theory, which allows real-time lattice studies and catches the quali¬ 
tative nature of the soft physics, which, has we have argued in Sec. 2, is of classical 
origin. In particular, they also examined the dependence on v of the Wilson loop 
(i.e. whether the approach to the light-cone from above or below is different for 
classical fields), finding it modest. 

Panero, Rummukainen and Schafer^^ have instead computed the soft contribu¬ 
tion to C{x_\_) and to q by means of lattice simulations of EQCD.° In the remainder 
of this Section, we will illustrate the basic principles of this treatment and its poten¬ 
tial application to other light-cone observables that are suited to it, such as q{6E), 
which will be defined in Sec. 5, or the asymptotic masses. To this end, we start by 
reviewing Electrostatic QCD (EQCD),®^^®® which is a dimensionally-reduced EFT 
of QCD describing the time-independent, long-wavelength A ~ l/{gT) modes of the 

° A different method, not based on the properties of Sec. 3, for estimating q on the lattice has 
been proposed. 
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latter theory. If treated perturbatively, it breaks down at very long wavelengths 
(very small momenta) A ^ l/{g^T) {k ~ g^T) where magnetic non-perturbative 
contributions®^’®® appear. If, on the other hand, it is treated non-perturbatively, 
both the electric and magnetic scales are treated correctly. Hence, a lattice calcula¬ 
tion in EQCD, such as the one of q we are considering, determines the soft {gT) and 
ultrasoft {g^T) contributions. Indeed, nonperturbative EQCD has been employed 
for many studies of the thermodynamics of QCD, as in this example,®® and of the 
gauge sector of the standard model, as in this review®^ and this recent work.®® 

The Lagrangian of EQCD reads 

C = + Tr ((AAo)®) + m|Tr (H®) + A 3 (Tr (H®))' , (72) 

where only operators of dimension up to three have been kept. Di = di — igE[Ai, •] 
is the standard covariant derivative, so that the Lagrangian describes a 3D SU(3) 
Yang-Mills field minimally coupled to an adjoint scalar field Aq, which has a mass 
and a quartic coupling (adjoint Higgs). Fermions, having only odd Matsubara 
frequencies proportional to ttT, are completely integrated out in the dimensional 
reduction. 

At leading order the matching between QCD and EQCD yields®^’®®’®® 

9e = g^{k)T, ml = g^{g)T^, A 3 = (73) 


where is the dimensionful coupling of EQCD and g{g) is the scale-dependent 
coupling of QCD. The mass parameter is equal to the leading-order Debye mass. 
The positive dimensionality of gE makes the theory super-renormalizeable. 

The EQCD action can be readily discretized for lattice studies We refer to^®’®® 
for details on this procedure. The implementation of the Wilson loop^®’®®’®® depicted 
in Fig. 3 is on the other hand somewhat non-trivial in the dimensionally-reduced 
theory, where the time dimension is integrated out. In the continuum, the Wilson 
line along the light-cone direction changes from 


u{xi 


■) = IP exp 


-ig f dx^ (A^{x~^) — A^{x~^)) 
J x'^ 


in Minkowskian QCD, 
(74) 


to 


U {Xz^f ^ ^z^i) 


-Pexp 


r^^.f 

-ig dxz{iAo{xz) - Az{xz)) 

9 ^z,i 


in EQCD, 


(75) 


where the A® in Eq. (74) is the temporal gauge field of Minkowskian 4D QCD, 
whereas the one in Eq. (75) is the adjoint scalar of Euclidean, 3D EQCD, the extra 
i coming from the Wick rotation. The other coordinates (xj_ and, in Eq. (74), x~) 
are understood to be constant and are not shown. 

The lattice version of Eq. (75) is also non-trivial, because in the discretized 
action the 3D gauge fields are naturally re-expressed as gauge links between the 
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lattice sites, while the Aq scalar field is defined on the lattice sites. This leads to, 
in the notation of^^ 

n^ —1 

L 3 (a;, aril) = ^ C /3 (x + anS) H [x a{n + 1)3) , (76) 

n—0 

where L 3 is the lattice equivalent of Eq. (75), with ni = Ija the length of the Wilson 
line in units of the lattice spacing a. C /3 is a gauge link in the 2 direction and i7(x) 
is a Hermitian, rather than unitary, matrix obtained by exponentiation of ^ 0 ( 2 ;), 

i.e. 

H{x) = exp[-a 5 |^o(a;)]- (77) 

Note that H (x) represents a parallel transporter along a real-time interval of length 
a. The transverse Wilson lines Li are instead the usual product of gauge links, so 
that, schematically, the Wilson loop takes the form depicted in Fig. 7. 







Fig. 7. The Wilson loop yielding C{x± = r) in EQCD. The dots on the lines in the horizontal 
direction represent the insertion of the Hermitian operator H at the lattice sites, whereas the lines 
connecting them are the gauge links in the z direction. Figure taken from the original reference.^® 

The lattice measurements^^ of C(x_l) through the above definition are summa¬ 
rized in Fig. 8 . As the authors remarked, once the value the Debye mass which 
includes the non-perturbative relative 0{g) correction®^ is used to set the scale and 
as input for the perturbative calculation of Caron-Huot, the agreement between lat¬ 
tice and perturbation theory becomes surprisingly good. If, on the other hand, the 
perturbative Debye mass is used, perturbation theory yields a significantly smaller 
C(x_l) at intermediate distances. 

The soft contribution to q can then be extracted from the short-distance behavior 
of C{x±) along the line of Eq. (37), i.e. ^^^ 4,64 

C(x_l) cix_L -k (?soft^ + 0{xl^). 


(78) 
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Coordinate-space collision kernel from EQCD 
(/i^=2, r=398 MeV, non-perturbative m^) 


Coordinate-space collision kernel from EQCD 
(n^= 2,1 ~ 2 GeV, non-perturbative 




Fig. 8. Results for C{xj_) at two different values of the temperature. The Debye mass used as 
a scale for the axes and plugged in the NLO perturbative calculation^^ is the non-perturbative 
one.®® Figure taken from the original reference.^® 

The qsoit extracted in this way also contains a residual logarithmical dependence on 
x±, reflecting the leading-order UV divergence of gsoft (see Eq. (66)) that is absorbed 
by the hard contribution. We refer to^^’®"^ for further details on the matching 
procedure and the addition of the leading-order hard contribution.^® Once this has 
been performed, Panero et al. report an estimate for g of 6 GeV^/fm at their lower 
temperature of 398 MeV, with an estimated uncertainty of 15 to 20%. We remark 
that higher-order contributions from the hard scale, as well as the possible collinear 
contribution mentioned in Sec. 3.3, are absent from this determination. 

A possible limitation to the approach we have outlined is that it relies on a 
separation of scale between the hard particles, with momenta of order ttT, and the 
soft and ultrasoft fields, characterized by gT and g^T. However, as the authors 
remark, the literature suggests (see for instance®®’®®’®® ) that analytical computa¬ 
tions relying on this separation of scales may be sufficiently accurate down to low 
temperatures, perhaps, suprisingly, down to®^ T ~ 2Tc. 

It is worth remarking that the NLO perturbative calculation predicts at the 
origin a negative linear slope , p i.e. Ci < 0 in Eq. (78), which is not observed 
in the lattice calculation. This can be attributed to discretization errors, which 
are more severe at short distances, corresponding to the UV region p± 3> mu 
in momentum space. At leading order in PT, the dominant UV behavior (\/p\) 
cancels between the longitudinal and transverse one-gluon exchanges, as shown 
in Eq. (65), leaving a m?jj/p\^ correction. This cancellation, while exact in the 
continuum, is only approximate on the lattice. D’Onofrio, Kurkela and Moore®® 
have estimated the associated error to be of order a/x±_ and hence especially relevant 
at short distances. Their computation of the renormalization properties of Eq. (76) 
to order a can help alleviating these discretization effects. 


PDue to the super-renormalizability of EQCD, each loop order causes a different power-law be¬ 
havior for C, so that higher orders cannot contaminate this effect. 
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We also note that, by comparing their results with the estimate of the magnetic 
contribution only to ^^33,64 can be performed similarly by simulating a Wil¬ 

son loop in the pure 3D gauge theory, Panero et al. have been able to establish 
that the magnetic contribution {g^T only) is subleading compared to the full non- 
perturbative EQCD determination, which includes the gT and g^T contributions. 

Future directions for this very promising new direction can include more re¬ 
fined measurements of 9 , for instance improving the approach to the continuum 
at short distances along the lines discussed before or extending the range of ex¬ 
plored temperatures. An altogether similar approach can be used to determine 
non-perturbatively the soft contribution to other operators that are amenable to a 
Euclidean treatment, such as q{6E) or the condensates appearing in the definition 
of the asymptotic masses (see Sec. 5.1). One could then envisage for the future 
the establishment of a factorization program for the computation of jet-related ob¬ 
servables: the soft physics, where pQCD struggles, is encoded into these effective 
operators measured on the lattice, whereas the harder scales are treated perturba- 
tively.'^ EFT techniques could be used to make the factorization well defined and 
deal with possible large logarithms. A hrst step in this direction has been taken in 
this study.It compares the screening masses of a specific correlator, extracted 
from lattice measurements, with their determination through a differential equation 
that is very similar to Eq. (43), up to a complex phase for C{x±). By using the 
lattice determination of the latter the authors find a very good agreement with the 
non-perturbative screening masses. 


5. Perturbative next-to-leading order 

As we have anticipated, the recent technological developments outlined in Sec. 3 
make the extension to next-to-leading order of the kinetic approach sketched in 
Sec. 2 possible. In that Section we have introduced the three fundamental processes 
that make up the leading-order collision operator, i.e. diffusion, large-angle scatter¬ 
ings and collinear splittings/joinings. Next-to-leading order contributions represent 
a relative 0{g) correction and come about because, as we have mentioned, soft glu¬ 
ons are highly occupied, with n^{q^ ^ gT) « T/q^ ^ !/(;. For this reason one can 
show that loops composed of soft propagators only are suppressed by a factor of g 
rather than 5 ^, giving rise to soft loop corrections. Another source of 0{g) correc¬ 
tions comes from regions of the LO calculation where a gluon becomes soft, without 
being properly treated, i.e. with HTL resummations. For instance, in large-angle 
scatterings, the integrations over the energies of the incoming and outgoing gluons 
stretch down to zero, including an 0{g) region of phase space where they are soft, 
but treated as massless, rather than HTL, excitations. At NLO one then needs to 
replace the improper LO evaluation of these mistreated regions with the proper one. 


^ A factorized approach sharing the same spirit, but using holographic methods for the description 
of the medium, has recently been presented in®^ . 
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In the following we will proceed to sketch the kinematical regions affected by 
NLO corrections and their evaluation. We will not dwell on mistreated regions, 
which, in a nutshell, require the identification and subtraction of the relevant coun¬ 
terterms in the corresponding points of the NLO calculations. 

Large-angle scatterings are not sensitive to soft loop corrections, as they would 
require internal and/or external soft lines, which are either excluded by construction 
or suppressed. Hard loops represent an 0{g^) correction, so that this region is not 
affected at NLO. 

Collinear processes are sensitive to soft loop corrections: they give rise to the 
0{g) correction^^ to C{x±) which we have mentioned in the previous Sections, as 
well as to an 0{g) correction^^ to the dispersion relations of the hard, collinear 
particles. 

Diffusion processes are also naturally sensitive to 0{g) corrections, since the 
typical momenta there are by construction soft. The soft loop corrections to 
have also been mentioned before, while those to will be introduced in the 
following. 

Finally, a new process appears at NLO. We label it semi-collinear. It is related 
to the example we made earlier on in this Section of a large-angle scattering with an 
incoming (or outgoing) soft gluon, which causes the angle between the two outgoing 
(or incoming) hard ones to be smaller than in the large-angle case but larger than in 
the pure collinear case. Indeed, as we shall see, this process will represent a bridge 
between the other three. 

Schematically, the collision operator then takes this form 

5C = JCcoll + iJCdiff [^t] -I- ^C'semi-collM, (79) 

where we anticipate, as noted by the g, dependence, that the last two processes will 
require regulation, but that the dependence on the regulator will vanish in the sum. 

5.1. The collinear region 

In the collinear regime the 0{g) corrections enter then in two places, as we antici¬ 
pated: both the effective thermal masses squared and the collision kernel C(a:_L) 
get 0{g) corrections''. 

For what concerns the thermal masses, an expression in terms of gauge-invariant 
fermionic and bosonic condensates supported on the light-cone is known.It reads® 

ml=g^CAZ, 

where v = (l,r>) is a null vector (|t;| = 1). Intuitively, Zg describes how thermal 
fluctuations of the gauge fields affect the propagation of a fast, light-like particle. 

For simplicity we do not illustate the two further mistreated regions'^® in the soft and semi- 
collinear limits respectively. 

®In QCD with quarks there is also a fermionic condensate. Both condensates contribute to the 
asymptotic masses of gluons and quarks. 
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It can be rewritten as 

— 1 

= dx+x+{v^Fr{x+)Uf{x+,Q)vXM^ ( 81 ) 

which makes the connection with q clearer: by comparing with Eq. (29), it is easy 
to see that they share the same Lorentz structure. At leading order Zg is dominated 
by hard {p ~ T) momenta and one easily has 





d^p 

{2-kYp 


«-b(p) 


T2 

T’ 


(82) 


The 0{g) contribution arises from the proper treatment of the IR limit of the above 
integration, requiring in principle HTL resummation. However, Zg is supported on 
a lightcone, so that the methods of Sec. 3 become applicable. As shown by Caron- 
Huot^° (further details can also be found here^^) the 0{g) contribution is easily 
obtained from the Euclidean zero-modes, yielding 


5Zg=T 


d^q 




{2t:Y {q^ + \q^ + m\g q 


-1 


1 


TruD 
2tt ' 


(83) 


so that the NLO correction to the asymptotic mass Sm^ reads 


Sm 


2 

OO 


9^Ca 


/ TniD 


(84) 


For what concerns the scattering kernel, Caron-Huot^^ has shown that three- 
body contributions vanish in the three-pole operator sketched below Eq. (43). 
Hence, also at NLO one just needs the sum of three two-body contributions at 
different transverse separations. The computation of these is nothing but the NLO 
contribution to C{x±). We have shown in detail in Sec. 3.3 how to “Euclideanize” 
the leading-order calculation. At NLO one then needs to compute one-loop dia¬ 
grams in three-dimensional EQCD, such as those shown in Fig. 9. The result for 



Fig. 9. Two sample diagrams necessary for the NLO evaluation of C{q±). The conventions are 
as in Fig. 4. 

SC{q±y^ and the corresponding 5C{x±_Y^ are not shown here for brevity. We point 
to the references for their expressions. We would like to stress again how the the 
reduction to the Euclidean, three-dimensional theory makes this calculation tremen¬ 
dously simpler. If we were to compute the diagrams of Fig. 9 in the HTL-resummed 
theory in four Minkowskian dimensions, we would have to deal with complicated 
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denominators of (up to 4) HTL propagators (see Eq. (54)) and the intricate HTL 
vertices. Indeed, in both our example diagrams we would need to consider them 
along the standard bare vertices. In the Euclideanized, three-dimensional theory, 
on the other hand, the propagators become simple massive or massless ones and the 
effective vertices just disappear, as they depend linearly on the external frequencies, 
which vanish in dimensional regularization.* 

In order to obtain the NLO correction to the collinear rate (40), Eq. (43) is 
then to be perturbed by and 5C. The hrst order in this perturbation 6f{h), 
is plugged in Eqs. (45) and (40). Details on this procedure may be found in these 
references.*"*’*® Algorithms for the solution of the differential equation with the 
correct boundary conditions are discussed in these other works. *^’*®’*’*’*"^ 

We furthermore stress that, should a factorized approach based on lattice mea¬ 
surements of soft observables, as outlined at the end of Sec. 4, be applied here, 
one would need a non-perturbative measurement of the aforementioned three-pole 
operator. It would be interesting to analyze the non-perturbative appearance of 
three-body forces in such an object. 

Finally, we remark that recently it has been pointed out*"'*’'’"* that collinear pro¬ 
cesses are sensitive to a large, double-logarithmic correction, which also shows up 
as a double-logarithmic correction to q. This term represent a relative 0{g^) cor¬ 
rection in our framework and would hence contribute to NNLO. It would arise from 
considering the first correction to the eikonal approximation in deriving Eq. (42). 


5.2. The diffusion sector at NLO 


As we anticipated, the Fokker-Planck diffusion equation (21) still applies. We then 
need the 0{g) corrections to q and q^, which unambiguously determine e at NLO 
through the procedure sketched in Eq. (17). Both can be calculated from their 
Wilson line definitions (26) and (27). The former is computed with the Euclidean 
technology of Sec. 3.2 and the latter with the sum rules of Sec. 3.3. However, 
before going in the details of the two separate calculations, one can proceed to a 
simplihcation,*®’^® which is attained by using the equation of motion of the Wilson 
line 0) = 0 to simplify to d^A~ (in all gauges but the A~ = 0 one) 

and by noting again that operator ordering is not relevant at NLO: classical soft 
gluons commute, as we have used in obtaining the simplified forms (29) and (30). 
The Wilson line structure then simplihes to 


p+oo 

J dx+{d-^A-‘^{x+)U'X\x+,0)d-^A-\0)), 

/ dx+{d+A-‘^[x+)Uf{x+,Q)d+A-\Q)). 

J — (x> 


(85) 

( 86 ) 


Their evaluation requires the computation of the diagrams shown in Fig. 10. 


^Indeed, there are no resummed vertices of the same order of the bare ones in EQCD, as we have 
shown in Sec. 4. 
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Fig. 10. Diagrams contributing to Sq and Sql at NLO. The blobs represent again the field strength 
insertions. Three- and four-gluon vertices are either bare or resummed HTL vertices. 


For what concerns q, the reduction to the Euclidean, three-dimensional the¬ 
ory introduced before makes this calculation too tremendously simpler. One then 
obtains for Sq, the NLO contribution to q^^ 

g^CnCAT^ mo 


Sq = 


(37r2-hl0-41n2). 


(87) 


327r2 T 

For what concerns Sq^, the sum-rule technology we have introduced in Sec. 3.3 
can be applied and leads to another tremendous computational simplification.^® As 
a general strategy, we remark that the Wilson line propagators depend only on the 
minus components of the momenta, so that the plus component, which we call q^, 
can be deformed again. This corresponds to expanding those diagrams for large, 
complex g"*". The leading contribution can be of order and the subleading 

one of order (g'*')”^. Higher-order terms are suppressed and can be neglected. 
The leading, 0((g+)®) term, once integrated along the contour, will give rise to a 
linear divergence, which has to match with the one arising from the corresponding 
mistreated region in the collinear process (see footnote r). Indeed, the calculation^® 
confirms that the two divergences cancel. 

For what concerns the 0{l/q~^) term, we already encountered such behavior 
at LO, where it gave rise to the asymptotic mass. In the photon case^^ it was 
found, rather surprisingly, that the NLO contribution to this term amounted to 
replacing the quark asymptotic mass with + Srn^, where Sm'^ is its soft 
0{g) correction, in the fermionic equivalent of Eq. (71) and then expanding to 
linear order in Sm^ (and hence in g). The explicit sum-rule computation^® of the 
diagrams in Fig. 10 yields the same (up to the different asymptotic mass) i.e. 




Smi 


mt 




Sml^ 




Sml 




ill 


\2 


0{g^), 


( 88 ) 


where Sm^ is given in Eq. (84). This result can be interpreted physically in the fol¬ 
lowing way: expanding for large (complex) g'*' takes the soft fields to approach their 
collinear limit, where the only effect of resummation is to introduce an asymptotic 
mass (at leading order) and the soft correction thereto (at next-to-leading order). 
For this same reason, HTL vertices become small and not relevant. 

^Caron-Huot obtained this result by applying Eq. (36) to his calculation of C{q±^ at NLO. The 
two methods are completely equivalent, as Eq. (26) can be formally derived from Eqs. (35) and 
(36) when C{q±^ is defined from the Wilson loop in Fig. 3, as shown in this reference. 
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Hence we obtain 


5qL = J 


cpqj_ qiSml^ _ g'^CnT6ml^ 
(2^)2(gi+m^)2 


47r 


In 




NLOt^ ' 


- 1 


(89) 


where we have introduced a regulator . As we will show, the semi-collinear 

region will remove the dependence on it, so that it should be taken to obey gT <C 
//Nlo << 


5.3. The semi-collinear region 

As we anticipated before, semi-collinear processes can be seen as 1 O 2 splitting 
processes where the opening angle (and hence the virtuality) are larger. Two exam¬ 
ples are drawn in Fig. 11. The scalings of this region are as follows: K ~ gT is soft. 


Q + K 





Q + K 



Fig. 11. Diagrams for two typical semi-collinear processes. In the first case the soft gluon is in 
the spacelike Landau cut, whereas in the second case it is on its timelike plasmon pole, represented 
by the black blob. 


whereas the two final-state particles are collinear, albeit with an increased virtuality 
and opening angle with respect to the collinear sector. The leading contribution 
then comes from < 7 + ^ T, g” ^ gT, ^ ~ gT^- 

Naive power-counting arguments would suggest that the semi-collinear region 
should contribute to leading order, as it is the largest slice of phase space where a soft 
gluon can attach to a 1 X-X 2 process. However, once all diagrams are summed and 
squared, a cancellation, first noticed in the context of photon radiation,^'' introduces 
an extra 0{g) suppression. Furthermore, since K ~ gT in all components, the 
contribution from timelike soft gluons, e.g. plasmons, is now allowed. This is 
contrasted by the collinear region, where kinematics enforce k~ ^ SE ~ g'^T 
k^,k±, thus restricting soft gluons to the space-like domain only. 

The contribution JCsemi-coU to the collision operator can be written in the same 
way as the collinear one, as given by Eq. (42), with the replacement of the collinear 
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rates with semi-collinear ones. In pure glue it reads 


16 


/ + 00 

doj 

-oo 


fp+<^ 


dT(j) + oj, uj) 


duj 


-fp0{P - 2w) 


semi—coll 

dT{p,u}) 


duj 


semi—coll 


(90) 


The derivation of the semi-collinear rates then requires the evaluation of processes 
of the form of Fig. 11, with p, g'*' q± ^ k±,k^. This results in the following 
factorized form^® 


dT{p,oj) 

dio 


semi—coll 


p'^Ca (1 - a: + 

np x{\ — x) 


(1 -I- n(w))(l -I- n{p 


uj)) 


where 


d?qx q{5E) 

( 2^)2 ql 

(91) 


5E = 


2u}{p — uj) ’ 


(92) 


and q{SE), which was first introduced for photon radiation,is a modified version 
of q which keeps track of the changes in the small light-cone components p~ and 
q ~, induced by the interaction with the soft degrees of freedom. 

In the field-theoretical language of Sec. 2, q{SE) is given by'' 


2^ poo 

q{6E) = ^ y dx+ {v^F^''^‘^{x+)Uf{x+,Q)vPFlM), (93) 


where we have made use of the simplification discussed in Eq. (85). The dependence 
on the minus components appears as a simple phase accumulation during the eikonal 
propagation, so that one can also easily see how Eq. (93) goes into Eq. (26) for 
5E 0. Eq. (93) can also be evaluated using Euclidean techniques, yielding^'^ 


q{5E) = g^CnT J 


d^k± 

(27r)2 [{kl + 6E^){kl + 6E^ + ml) 




2SE^ 

kl+6E^ 


(94) 


As we have anticipated at the beginning of this Section, a counterterm from the 
large-angle scattering region needs to be subtracted, as well as a collinear one (see 
footnote r) . The former arises from taking the mo —t 0 limit of Eq. (94) (there is 
no resummation taking place on external legs at LO) and the latter by taking the 
SE —>■ 0 limit (changes in the minus component are not accounted for in the collinear 


'^Strictly speaking, this is the version that applies to photon radiation and in this case we would 
need a more complicated three-pole expression, along the lines of the discussion that followed 
Eq. (43). However, we are evaluating it at the leading-order, one-gluon exchange level, so that the 
more complicated form is not necessary. 


























February 13, 2015 1:25 


World Scientific Review Volume - 9.75in x 6.5in 


eloss 


Parton energy loss and pT broadening at NLO in high temperature QCD 33 


limit). Once these subtractions are performed, the semi-collinear rate becomes 

dVip^co) g^ClT^^ , , , . ^^ /■ f 1 

-( 1 + n(a;))(l+ n(p-a;)) J J 


dio 


semi—coll 


irp 


( 2^)2 ql 


(1 - X + 




:(1 — x) 


rnjjk\ 


‘■D 


{kl+SE^){kl+dE^+ ml,) 


kl 


‘‘D 


(95) 


We stress that the semi-collinear collision operator is given by Eq. (90). 

The qx integration in Eq. (95) is to be understood as IR-regulated by . 
One can show^® how the small-a;-and-g_L region gives rise to IR logarithms that 
cancel the dependence of diffusion processes. Details of how the transverse 

integrations can be carried out analytically are also given there. The w integration 
remains to be performed numerically. 


6. Summary and conclusions 

In this review we have shown how at leading and next-to-leading order the propa¬ 
gation, energy loss, and momentum diffusion of high-energy particles in the QGP 
can be cast in the form of a Boltzmann equation describing the interactions be¬ 
tween the jet particles and the hard and soft constituents of the plasma (Eq. (1) 
and Eq. (79)). In Sec. 2 we introduced this formalism at leading order, where three 
processes contribute. These processes include 2 e-)- 2 scatterings of the jet with the 
hard particles, drag and (longitudinal and transverse) momentum diffusion induced 
by the soft highly occupied gluonic background, and collinear bremsstrahlung (and 
merging) induced by the transverse momentum diffusion. The 2 o 2 scatterings 
are dealt with in a standard kinetic fashion (Eq. (10)), the collinear splittings and 
joinings give rise to a somewhat simpler rate term (Eq. (42)), and the interactions 
with the soft background are described by a Fokker-Planck equation (Eq. (21)) or 
equivalently by a Langevin equation (Eq. (23)). 

We have furthermore shown how the contributions of the soft degrees of freedom 
to the kinetic picture (the diffusion coefficients in the Fokker-Planck equation and 
the transverse scattering kernel in the collinear rate) can be cleanly defined as 
correlators of Wilson lines, or Wilson loops, supported on the lightcone or on a 
light front (see Eqs. (26), (27) and (39)). These Wilson lines are to be evaluated 
using the Hard Thermal Loop effective theory. We devote Sec. 3 to the illustration of 
the recent understanding of thermal field theory (including Hard Thermal Loops) 
on the lightcone. We show how the causal properties of amplitudes at light-like 
separations bring about a dramatic simplification in their calculation.^^ To put 
it simply, jet and hard particles, which move at ~ c, probe the medium in a 
very simple (thermodynamical) way, since the soft background “does not have the 
time” to respond to the perturbation they cause. At a practical level, we show 
how the computation of certain lightlike correlators, such as transverse momentum 
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diffusion, “euclideanizes”, that is, it reduces to a much simpler calculation in the 
three-dimensional Euclidean theory EQCD (see Eqs. (64) and (65)). We also show 
how a different set of lightlike correlators, such as longitudinal diffusion, does not 
euclideanize. This case also simplifies, however, and the same causal properties on 
the lightcone make it sensitive only to the thermal mass and soft corrections thereto 
(see Eqs. (70) and (71)). 

In Sec. 4 and Sec. 5 we reviewed two major consequences of these lightcone 
simplifications. First, we showed how operators that euclideanize can be measured 
on the lattice. Second, we reviewed the consequences on perturbation theory. At 
leading order the lightcone simplifications result in simple analytical closed forms 
in place of relatively straightforward numerical integrals. At NLO the impact is 
much more far-reaching - calculations that would have been extremely intricate, if 
not impossible, can now be performed with relative ease. 

In more detail, in Sec. 4 we reviewed the basics of EQCD, and described the first 
non-perturbative calculations of q and the transverse momentum scattering kernel 
using this effective theory. We would like to emphasize how these first results 
open a new avenue of research - all other Euclidean operators can be measured on 
the lattice in the same way. This creates the tantalizing possibility of a factorized 
approach to kinetics, where perturbation theory is used at the hard scale to compute 
the 2 o 2 scatterings and the splitting processes, while the 3D lattice is used at 
the soft and ultrasoft scales to compute the diffusion processes and the transverse 
scatterings leading to bremsstrahlung. 

In Sec. 5 we give an overview of the NLO generalization of the LO Boltzmann 
equation. We showed how the three processes at LO are corrected at NLO (see 
Eq. (79)). Besides the corrections to the LO diffusion and collinear processes, a 
new semi-collinear process enters, which appears as a bridge between the others, 
being collinear but with relaxed constraints, i.e. going beyond strict collinearity. 
In the evaluation of each of these processes we sketched how the lightcone sim¬ 
plifications described above make the NLO analysis possible. We described the 
NLO calculations of transverse^^ and longitudinal^® momentum diffusion (Eq. (87) 
and Eq. (89)), and introduced two related Euclidean operators: the condensate Zg 
responsible for the thermal mass (Eq. (80)), and q{dE), a modified version of q 
relevant for the semicollinear bremsstrahlung (Eq. (93)). The corrections to the 
longitudinal diffusion and drag coefficients at NLO are due to corrections to the 
thermal mass (or Zg). Finally, we presented complete expressions for the NLO 
rates in the Boltzmann equation for pure gauge theory. 

The extension of the NLO Boltzmann equation to include light quarks requires 
one extra process at leading- and next-to-leading order, the conversion process, 
which can be seen as a fermionic analogue of diffusion.^® The interaction of hard 
or jet particles with soft quarks results in a conversion from a quark to a gluon 
with approximately the same momentum and vice-versa. The computation of these 
rates involve fermionic correlators on the light cone. It turns out that the light- 
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cone simplifications of the fermionic HTLs are strikingly similar to the longitudinal 
diffusion described in Sec. 3.3, and consequently the conversion rates are sensitive 
only to the (fermionic) asymptotic mass and its soft corrections. These corrections 
to the fermion mass are computable with the 3D effective theory, in much the way 
that the effective theory can be used to compute corrections to Zg. 

One very important question that we leave to future works is the impact of these 
developments on calculations of jet energy loss and their comparison to experimental 
data. We remark that the Monte Carlo event generator MARTINF^ currently im¬ 
plements a kinetic approach that is very close to the leading-order picture sketched 
in Sec. 2. It is then an ideal candidate for the inclusion of the developments discussed 
in Sec. 4 and Sec. 5. The reorganization in terms of 2 y-)- 2 scattering, diffusion and 
collinear processes, as well as the implementation of the NLO corrections discussed 
in Sec. 5, is already underway. This can be easily complemented by the inclusion 
of non-perturbative input. Besides the measured transverse scattering kernel, lat¬ 
tice calculations for q{5E) and could easily be included into MARTINI, should 
they become available in the future. It would also be interesting to see the results 
of this numerical implementation for the angular structure of jets, for which recent 
order-of-magnitude estimates from perturbation theory are available.^® 

As a general remark, we find it difficult to tell a priori the impact the NLO 
corrections of Sec. 5 will have. The recent calculation to NLO of the thermal 
photon rate,^^ (which includes many of the ingredients presented here, such as Eu- 
clideanizations, light-cone causality, collinear and semi-collinear processes) showed 
how the NLO corrections naturally divide into two large, and largely canceling, con¬ 
tributions - see Fig. 12. Clearly the NLO corrections are modest over a significant 
range of photon momentum. The positive contribution arose from collinear pro¬ 
cesses enhanced by a NLO increase of the scattering kernel, and a NLO decrease of 
the thermal mass, while the negative contribution arose from the semi-collinear and 
NLO soft processes. The large cancellation between these two corrections is largely 
accidental and is dependent on the details of the medium, such as the number of 
colors and flavors. Nevertheless, we anticipate a similar cancellation for the current 
more complicated case of parton energy loss. 

Finally, we note that the kinetic approach we have outlined is often used to 
compute the transport coefficients of QCD plasmas.While an extension to NLO 
of these calculations requires more than what is outlined in Sec. 5, a computation 
of transport coefficients in thermal QCD beyond leading order will almost certainly 
make considerable use of the euclidean light-cone simplifications described in this 
review. 
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Fig. 12. The differential photon emission rate dV^/dk relative to the leading order rate as a 
function of the photon momentum k/T at a fixed coupling Qs = 0.3 in a Nf = 3 QGP. The 
full next-to-leading order rate^^ (LO+NLO) is a sum of the leading order rate (LO), a collinear 
correction (coll), and a soft-Fsemi-collinear correction (soft-Fsc). The dashed curve labeled LO-Fcoll 
shows the ratio of rates when only the collinear correction is included, with the analogous notation 
for the LO + soft-Fsc curve. Figure taken from the original reference. 
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